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Star formation is regulated through a variety of feedback processes. In this study, we treat feedback by X-rays and discuss its impli- 
cations. Our aim is to investigate whether star formation is significantly affected when a star forming cloud resides in the vicinity of 
a strong X-ray source. We perform an Eulerian grid simulation with embedded Lagrangian sink particles of a collapsing molecular 
cloud near a massive, 10 7 M Q black hole. The chemical and thermal changes caused by radiation are incorporated into the FLASH 
code. When there is strong X-ray feedback the star forming cloud fragments into larger clumps whereby fewer but more massive 
protostellar cores are formed. Competitive accretion has a strong impact on the mass function and a near-flat, non-Salpeter IMF 
results. 
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1. Introduction 



At a distance of 100 pc, a massive 10 7 M bl ack hole is able 
to em it X-rays with a flux of 100 erg s _1 crrT 2 dMeijerink et al.l 



2007). The presence of X-rays in star forming regions alters fun- 
damental processes in star forming molecular clouds, that even- 
tually determine stellar masses, like ionization and gas and dust 
heating, and do s o up to column densiti es of 10 24 cm~ 2 and dis- 



l_L I tances of 300 pc dSchleicher et al.ll2010h . 
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In our local neighborhood, the initial mass function (IMF) 
is observed to be a power-law, nicely following a Salpeter 
slope. The IMF is assumed to have a universa l shape and 
spatia l ly well resolved stud ies dChabried 120031: ISabbi et ail 
l2007t lElmegreen et all l2008h confirm this at low and high 
masses. When we turn to observations of more radical en- 
vironme nts like the G al axy center, hints for different IMFs 
appear (iFigerl I2005allbt iPaumard et .al. 20061: lEspinoza et al.l 
l2009t lElmegreenll2009t iBartko et al.ll2010t). However, the evi - 
dence is also claimed to be controversial ( Bastian et al.l l2010b . 
Nuclei of active galaxies, e.g., ULIRGs l ike Arp 220 and 
Markarian 231, see van der Werf et al. 1 (12011 . enjoy conditions 
much more extreme than the Milky Way galactic center. Little 
is known on the IMF there, so theory and simulations are 
needed to guide understanding. The idea arose that massive 
stars form more easily i n these extreme regions and that the 
IMF might be top-heavy (|Brassington et al ] l2007tlKlessen et al.l 



120071: iDabringhausen et al. 1 120091) . Other explanations are possi- 
ble, still, it is worthwile to test this hypothesis numerically. 

Until now, there have been many studies of feedback pro 
cesses and external influences on star formation {Bonnell & Ric 



2008; Wada 2008; Wada et al. 2009r. iKrumholz et aLl^OlOt lBafr 



2010). However, no numerical simulations have been performed 
yet of a collapsing molecular cloud under strong incident X-ray 
radiation. In this paper, we present the first 3D numerical simu- 
lations on the effects of X-rays on star forming molecular clouds 
that reside close to a massive black hole. In the next sections, we 
compare such a molecular cloud to an X-ray free environment 



and evaluate the star formation efficiency, stellar masses, and the 
resulting IMF. 

2. Numerical model 

2.1. The FLASH code 

All simulations in this study have been performed with 
the Eulerian hydrodynamical and N-body code FLASH 



dFrvxell et al.ll2000t iDubev et alJl2009l) . We used FLAS H with 
the im proved gravity solver that was added in version 3 dRickeii 
I2008I) . FLASH is a modular based, strongly-scaled parallel code 
that is specialized in adaptively refined meshes in which one can 
use particles in conjunction with the grid. 

The hydrodynamic e quations are solved using the piecewise 



parabolic method (PPM, IColella & Woodward 1 1984. which i s 



an improved version of Godunov's method (lGodunov|[l959l) . 
PPM is particularly well suited for flows involving disconti- 
nuities, such as shocks, that are strongly present in this study. 
FLASH is provided with many and extensively tested modules 
that encompass a broad range of physics. Many of these mod- 
ules, including hydrodynamics, thermodynamics, (self-)gravity, 
particles, turbulence, and shocks are standard ingredients for in- 
terstellar physics and star formation and are thus incorporated 
in this work. Several additions were made to the code to follow 
the (radiation)physics in detail, such as sink particles, radiative 
transfer, multi-scale turbulence, and refinement criteria based on 
Jeans length and sink particles. The main additions are described 
in the following subsections. 

2.2. Sink particles 

Sink particles are point particles that can grow in mass by accret- 
ing gas and merge with other particles, however, they cannot lose 
mass or fragment. These particles represent compact objects, in 
our case protostars. Sink particles are a necessary ingredient if 
one wants to follow a density evolution over a large dynamic 
range, like a collapsing cloud. We created a sink particle module 
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that takes care of the particle creation, determines their masses, 
handles Bondi-Hoyle type accretion, and tracks the particles if 
they are eligible for merging. 

A sink particle is created when it passes several criteria for 
indefinite collapse as given bv lFederrathetalJd2010l) . The mass 
it obtains at creation is determined b y a density threshold th at is 
set by the maximum resolution. The iTruelove et al.l d 1997b cri- 
terion states that the Jeans length should not be resolved by less 
then 4 cells if one wants to avoid artificial fragmentation. This 
criterion can be rewritten in the form of a density limit. Any ex- 
cess density is taken away from the grid and added to the mass 
of the particle. The sink particl e algorithm created for this study 
follows the method created by iRrumholz et al. (2004). For any 
details of this method, we refer the reader to that paper. 



2.3. Radiative transfer 

In order to update the effects caused by X-rays on the temper- 
ature of the gas, we ported an X-r ay dominated region chem- 
ical code (XDR code) created by iMeiierink & Spaans] d2005l) 
into FLASH. This code incorporates all of the heating (photo- 
ionization, yielding non-thermal electrons) and cooling pro- 
cesses from atomic (fine-structure, semi-forbidden) and molec- 
ular transitions (CO, H2, H2O). Effects from internal UV, cos- 
mic rays, and dust-gas coupling are treated as well. Given an 
X-ray flux, gas density, and column density along the line of 
sight to the source, the XDR code calculates the temperature and 
the chemical abundances. This output is fed into the simulation 
at every iteration. Most of the computation is spent finding the 
column densities for every cell. A ray-tracing algorithm, specif- 
ically created for this purpose, searches the grid and sums up 
the column densities of each cell lying along the line of sight 
from the source, the accreting black hole, to the target cell. The 
X-ray flux is an E~ 09 power law between 1 and 100 keV. X-ray 
scattering is not very important, but is nonetheless treated in the 
XDR-code. A uniform background of cosmic rays prevents the 
temperature from dropping below 10 K. For this, a cosmic ray 
ionizatio n rate typical for the Milky Way £ = 5 x 10~ 17 s , is 
assumed (ISpaans & Silkll2005l) 



2.4. Initial conditions 

We create a spherical gas cloud with solar abundances at a dis- 
tance of 10 parsec from a 1O 7 M black hole. We run two sepa- 
rate simulations and name them simulation A and B. Simulation 
A represents a molecular cloud near an active black hole under 
the impact of X-rays. Simulation B has a cloud near an inactive 
black hole and has isothermal conditions, with an equation of 
state of the form P oc p. The XDR code that updates the thermo- 
dynamics is only linked with simulation A where we do have X- 
rays. The 1O 7 M black hole, at 10% of Eddington, yields a flu x 
of 160 erg s" 1 cm' 2 with some extinction (IMeiierink et al.l2007l) . 
IWada etal.l(l2009l) show that column densities of 10 225 cm 2 can 



5.0e+24 



typically exist in the central R 10 pc of an AGN, leading to a 1 
keV optical depth of 2-3. Furthermore, Wa da et all ((2009) show 
that column densities as large as 10 24 cm~ 2 occur and persist in a 
statistical sense in the dynamically active inner 20 pc. A gaseous 
cloud can thus be shielded by this clumpy medium around the 
AGN and remain cold, while the temperature rises rapidly once 
the cloud is exposed to the radiation. Both our simulations start 
with the same initial conditions, but we expose simulation A to 
an X-ray source once the simulation starts. In this, heating is 
nearly immediate, since the timescale for heating is much shorter 
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Fig. 1. Column density (cm -2 ) slice of simulation A at t = tff . 
The colors represent the column density along the line of sight 
to the black hole, which is located at the upper left side (the 
black arrows indicate the direction of radiation). The density is 
shown in black contours. Contour levels are 1, 4, 16, 64, 256, 
1024 xlO 4 cm 3 . The white spheres display the location of the 
sink particles and the red arrow indicates the direction of the 
cloud's orbital motion. 



than the collapse time, theat ^ tff, and of the order of 10 _1 years. 
All other conditions are the same for both simulations. 

The simulations are set up with an initial random, 
divergence-free turbulent velocity field and a characteristic 
FWHM of 5 km/s that agrees well with molecular clouds ob- 
served in active regions dPerez-Beaupuits et al.ll2009l) . These are 
supersonic flows with Mach numbers of up to 25, where the 
isothermal sound speed of the cloud is c s = 0.19 km/s (for 
T=10K) and can go up to a maximum of 5 km/s when the cloud 
is heated by X-rays (T^IO 4 K). We do not drive the turbulence 
but follow its decay. The turbulence is applied over all scales 
with a power spectrum of P(k ) oc k~ 4 following the empirical 
laws f or compressible fluid s ( Larson|[l98U iMvers & Gammid 
11999b iHever & Brunl l2004h . We start the simulations with a 
cloud that is in a stable Keplerian orbit around the black hole. 
Shear that is introduced by the black hole is taken into account. 
The maximum velocity difference imposed by the black hole, 
Av s heai = 2.2 km/s across the cloud, is of the order of the ap- 
plied initial turbulence. This process keeps the turbulence strong 
to large dynamical times, i.e 



> 1 t ff , with tff = yj3n/32Gp 
and 10 5 years3 throughout this work. The shearing time, t s h ea r = 
Dcioud/Av shear, is almost 3 times larger than the cloud free-fall 
time and gravitationally bound (roughly) spherical clouds are 
likely to exist at densities of ~10 5 crrr 3 . 

The molecular cloud starts with a uniform number density 
of 10 5 cm~ 3 and has a size of 0.33 pc in radius. With a mean 
molecular weight of p. — 2.3, the total mass of the cloud amounts 
to 800 solar masses. The rest of the medium is filled with gas that 
has a uniform density of 100 crrr 3 . The simulation box, a cube 
of size 24 pc, has outflow boundaries and is isolated in terms 
of gravity. We increase the resolution where needed according 
to a self-developed Jeans criterion. The algorithm calculates the 
Jeans length at every grid cell, compares it against the Truelove 
criterion, and adds resolution when this is about to be violated. 
The maximum grid resolution that we allow for any simulation is 



Given the initial condition of p = 3.84 • 10 l9 g/cm 3 
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Fig. 2. Density morphology of simulation A at t=0.4tjf in 2D slices through the center of the cloud with axes in parsec. The color 
represents the number density (cm -3 ). On the left, an XY slice is shown with arrows representing the direction of radiation emanating 
from the black hole, which is located at the upper left side. On the right, an XZ slice is shown with radiation that is coming from 
the left and under an angle of 16 degrees with the slice. 



8 192 3 cells. With the box size of 24 pc, the maximum spatial res- 
olution becomes 8.8 X 10 15 cm. If gas continues to collapse and 
needs higher resolution beyond the maximum refinement level, 
sink particles are created thereby taking density away from the 
gas such that mass and momentum is conserved and the resolu- 
tion criterion is not violated. We give sink particles an accretion 
radius of effectively 3 cells, that is, 2.6xl0 16 cm 1760 AU. The 
temperature in both simulations is initialized at 10 K, but for the 
simulation with X-rays, the temperature is updated by the XDR 
code, and changes quickly (<10~ 4 tff = 1 timestep) after initial- 
ization. 

3. Results 

3.1. X-rays versus no X-rays: IMF and Jeans mass 

Both simulations, A and B, are followed for 2xl0 5 years. This is 
approximately two free-fall times for a 10 5 cirT 3 cloud. X-rays 
can heat the cloud to as high as 6000 K at low column densities 
(< 10 21 cm -2 ), this is the case for the side that is directly exposed 
to the X-ray source, but also to as low as 10 K at high column 
densities (> 10 24 cm -2 ). Fig.[T]shows a column density plot to 
this effect. 

The directional heating increases the pressure and causes 
the gas to expand and evaporate on the irradiated face of the 
cloud. The molecular cloud is compressed, loses mass, and an 
ionizing pressure flow travels inward. We see that this compres- 
sion creates a density increase of about half an order of mag- 
nitude within a free-fall time as compared to simulation B. The 
conical compression front is disrupted where the turbulence cre- 
ates sub-pc scale gaps (0.01-0.05 pc) and radiation is able to pen- 
etrate, as is illustrated in Fig. [2] Those regions are also heated up 
and pressurized. This causes the pressure front of the irradiated 
side of the cloud to break up. We see finger-like shapes forming, 
with a high density head, and the gas that is lying in its shadow is 
well shielded and very cold (about 10 K). The increased density 
induces star formation. We find that sink particles are created in 
the compressed cloud edge much earlier than in the shielded and 



colder parts. A phase diagram is plotted in Fig.[3]that shows the 
decline of temperature as the density increases. The secondary 
band with a steep decrease in temperature at low densities is 
the direct result of shielding. We note that the reason that fewer 
points appear in the plot in the shielded regions is merely due to 
the difference in resolution determined by the refinement crite- 
ria. 
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Fig. 3. Temperature-density phase diagram of the X-ray irradi- 
ated simulation (A) after 1 dynamical evolution, t = tff. 



Simulation A starts forming protostars at around two-thirds 
of the initial free-fall time, which is later than simulation B, but 
at about three-quarters of the free-fall time, we see a sudden in- 
crease of star formation. Simulation B, on the other hand, has a 
gradual increase of protostars, starting at about a half tff, only to 
peak close to the free-fall time of 10 5 years. In the end, simula- 
tion B has created more sink particles, 153 against 118, but their 
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Fig. 4. The initial mass function of the simulations at t = tff. Left panel shows the IMF for simulation A, and the right panel shows 
the IMF for simulation B. The green dashed line represents the Salpeter function. The blue dot-dashed line represents the Chabrier 
IMF, as fitted to simulation B. A power-law fit is applied to the data above the turn-over mass, ~1.2 M left and ~0.7 M right, and 
is shown as the solid purple line. 



masses are much lower. We plot the particle masses against the 
number of particles in Fig. [4] 

The plots show that the X-ray dominated environment 
clearly has a higher characteristic mass as well as a higher min- 
imum and maximum mass compared to simulation B. Higher 
sink particle masses are seen consistently throughout simulation 
A. Another difference is evident when looking at the slope of the 
IMF, which is defined as 



dN 
dM 



oc M 



dlogN 
dlogM 



-a + 1 = r, 



(1) 



with N the number of stars in a range of mass dM, a the power- 
law index, and F the slope above the characteristic mass of ~0.3- 
0.5 M . 

Simulation A has a much flatter slope than Simulation B, 
which is Fa = -0.28 after 10 5 years of evolution, flatter than 
the Salpeter slope of -1.35 dSalpeterlll955l) . There is also a steep 
decrease in the number of stars above ~3 M . This slope is very 
steep but uncertain and ca n depend strongly on late accretion and 
merging dDib et alJl20loh . Simulation B, has a slope of Tb = - 
1.37 which, in fact, is very similar to the Salpeter slope. Our 
simulations show that these trends for simulation B are main- 
tained even at higher dynamical times, t = 2tff, in agreement 
wi th the time i ndepe ndent, competitive accretion driven results 
of lClarketal1(l2008l) for initial energies of |E grav | > Eyn- We at- 
tribute the difference in F partly to the higher gas temperatures 
(~50 K) caused by X-rays, which in turn increases the Jeans 
mass and the fragmentation mass scale. Most of the massive sink 
particles are created in the dense fragments and trace the warm 
X-ray heated (T>50 K) molecular gas. The Jeans mass in case B 
is about 0.2 M after one free-fall time, whereas the Jeans mass 
in case A can be much higher, but typically ranges from 0.5 to 
2 M in the > 10 6 cm 3 gas. T he Jeans mass is approximately 
given by dFrieswiik et al.ll2007h ; 



M, = U 



i 90 3 i 
p -2 ^—Ji n-s (Mo). 



(2) 



In response to the high temperature in the XDR gas, only 
high density regions collapse creating more massive protostars 



and, in the end, fewer of them. Furthermore, the cloud is shielded 
from radiation with increasing column density, causing less X- 
ray penetration and thus less heating, lowering the Jeans mass. 
As a result, gas becomes more compressible and the density in- 
creases, starting a snowball effect. This causes a steep decrease 
of the Jeans mass, giving rise to a sharp increase in star forma- 
tion. Another feature of the aforementioned effect is that almost 
all of the particles that are created in simulation A, are formed 
in the high density front, the finger head, of the cloud, where the 
snowball effect is initiated. 



3.2. Sink particle behavior 

We let the code check whether sink particles come too close to 
one another and if they should merge. This is the case if the ve- 
locity difference between two particles is less than the escape ve- 
locity and if the merging time is shorter than the simulation time 
step. Merging is not often seen. We define the merging time as 
half of a degenerate orbit time, t merge = (8G(mi +m2))~^ 2 , 
where mi and iri2 are the masses of two independent sink parti- 
cles and Rn is the distance between the two. 

Mergers can potentially increase the mass considerably, 
but most of the protostellar mass gain comes from accre- 
tion. To this end, we a dopted a Bondi-Hoyle type of accretion 
(iRrumholz et al.l 120041) . This type of accretion arises when a 
homogeneous flow of matter at infinity moves non-radially to- 
wards the accretor. Accretion increases with protostellar mass, 
but drops with decreasing ambient density and higher Mach 
numbers. We see that the ratio of the most massive to the sec- 
ond most massive sink particle in a local cloud region grows 
in time. With this, our expectation is confirmed that protostars 
compete with each other for material and that high mass stars 
as well as stars that lie in deep pot ential wells accrete more 
dBonnell & Batdl2006tlBonneH2008l) . 

The total mass of the particles after one free-fall time for 
simulation B is 104 M , i.e., about l/8th of the initial cloud 
mass. Quite interestingly, simulation A is able to convert more 
gas mass into stars, 230 M . This is a strong argument for ef- 
ficient X-ray induced star formation in AGN. The star forma- 



S. Hocuk & M. Spaans: The impact of X-rays on molecular cloud fragmentation 



5 



tion efficiencies mentioned here are upper limits, since these 
calculations do not contain feedback effects from young stars. 
Feedba ck effects, like ou tflows, should decrease the efficiency 
overall dWang et alJIioToh . 

4. Conclusions and Discussion 

We have performed two 3D simulations of similar molecu- 
lar clouds, each at 10 pc distance from a supermassive black 
hole and followed their evolution. In one case (A), we ex- 
pose the cloud to an active black hole, producing a strong, 160 
erg s crrT 2 , X-ray flux. In the other case (B) the black hole was 
inactive and the molecular cloud had isothermal (10 K) condi- 
tions throughout the run. We saw clear differences between the 
simulations. 

For the X-rays included run, we find that the molecular 
cloud is heated at the irradiated side and an ionizing pressure 
front is formed. This conic pressure front breaks up due to tur- 
bulent motions, forming finger-like structures that can be seen 
in column density plots, see Fig.Q] The density increases at the 
top of the compression front and the shielded gas cools. Star for- 
mation is initially delayed due to the higher temperatures, but 
the cloud continues to contract and after the critical density is 
reached sink particles start to from. Since the temperature in the 
shielded parts of the cloud decreases with increasing density, the 
Jeans mass drops at an amplified rate. Consequently, protostar 
formation increases sharply around 0.75 tjf. Despite this, fewer 
sink particles are created in total with respect to simulation B, 
but they are more massive. The latter is a consequence of the 
high temperature (~50 K) and Jeans mass (>1 M ) in X-ray ir- 
radiated gas. In the end, case A ends up with more total stellar 
mass after one free-fall time. These protostars also accrete more 
material as the accretion rate scales with mass and density. This 
becomes increasingly more important for the mass growth due to 
the deeper potential well created and favored by massive stars. 
The two effects together cause that competitive accretion domi- 
nates the mass growth and strongly affects the shape of the mass 
function. The resulting IMF has a higher characteristic mass and 
a near-flat, non-Salpeter slope. We summarize the main points 
by stating how case A behaves with respect to B: 

- Protostars are created at a later stage, around 0.65 tff (0.45 tg 
for B). 

- Fewer protostars are created but they have higher masses, 
M=0.8-6 M (0.2-3.6 M Q for B). 

- Although protostars are created later, there is a burst mode 
around 3/4th of the collapse time of 10 s years. 

- The total stellar mass formed from the gas is more than twice 
higher after 10 5 years and the efficiency is about 28% (13% 
for B). 

- Competitive accretion is more influential in shaping the IMF. 

- The IMF has a much flatter slope, T = -0.28 (-1.37 for B). 

- The characteristic mass of the IMF is higher, M c h ar ~1.2 M Q 
(-0.7 M Q for B). 

All of these results tell us two main things. Firstly, star for- 
mation is induced by X-rays and massive stars can form in this 
way with high efficiency. Secondly, the resulting IMF differs 
from a Salpeter shape. 

The pressure force exerted by radiation from the black hole 
is found to be modest. The bulk of the radiation pressure will re- 
sults from the UV, since it dominates the bolometric luminosity 
L for very massive black holes. Typically, the change in momen- 
tum of the gas is proportional to tfir L/c, for the speed of light 



c and mean FIR optical depth tfir. The latter comes in because 
the dust, if it absorbs the bulk of the hard radiation, re-emits this 
in the FIR (100-300 fim). The ratio of radiation pressure over 
thermal pressure, for a typical tfir of ~0.1, is about 10% at 10 
pc from the black hole 

We did not include any model for stellar feedback in this 
study. The total time of the simulations is not long enough for 
stellar evolution effects to play a major role, however, outflows 
from young stars can produce significant winds that disrupt gas 
build-up. This may affect the later time (> 1 tjf) accretion onto 
protostars. 

There is also strong gravitational shear, which helps to drive 
the turbulence but also compresses some regions in the center 
of the cloud along the direction perpendicular t o the orbit. This 
increa ses the star formation efficiency overall dBonnell & Ricd 
I2008h . but does not affect the conclusions when comparing the 
two cases with each other. 

There are several other parameters that can be studied for 
their effectiveness, like the X-ray flux, turbulent strength, cloud 
size and initial gas density, distance to, and mass of, the black 
hole. In a follow-up paper we plan to present a detailed param- 
eter study to this effect including the impact of UV and cosmic 
rays. 
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